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Am J Physiol Lung Cell Mol Physiol 280: L841-L852, 
2001. — Epithelial cells lining the airways are thought to play 
a prominent role in respiratory diseases. We utilized cDNA 
representational difference analysis to identify the genes in 
which expression is induced by the proinflammatory cyto- 
kines tumor necrosis factor-a and interleukin-ip in primary 
human bronchial epithelial cells and hence are relevant to 
airway inflammation. Hybridization of the subtraction prod- 
uct to arrayed cDNAs indicated that known tumor necrosis 
factor-a- and interleukin-ip-inducible genes such as B94, 
Zfp36, and regulated on activation normal T cell expressed 
and secreted were represented, confirming the success of the 
subtraction experiment. A 1,152-clone library potentially 
representing genes with higher transcript levels in cytokine- 
treated human bronchial epithelial cells was generated and 
sequenced. Sequence similarity searches indicated that these 
clones represented 57 genes of known function, 1 gene of 
unknown function, 6 expressed sequence tags, and 2 novel 
sequences. The expression of 19 of these clones was studied 
by a combination of Northern blotting and RT-PCR analyses 
and confirmation of differential expression for 10 known 
genes, 2 expressed sequence tags, and a novel sequence not 
represented in any of the public databases was obtained. 
Thus cDNA representational difference analysis was utilized 
to isolate known and novel differentially expressed genes, 
which putatively play a role in airway inflammation. 

representational difference analysis; human bronchial epi- 
thelial cells; tumor necrosis factor-a; interleukin-l(3; airway 
disease 



the airway epithelium has an important role in many of 
the responses seen in respiratory diseases. It consti- 
tutes the interface between the internal milieu and the 
external environment. It is the first point of contact for 
inhaled substances, respiratory viruses, and airborne 
allergens. It is also a primary target for inhaled drugs. 
The airway epithelial cells are classically considered to 
provide the physical barrier between the external and 
internal environments. A recent study (1) has indi- 
cated that airway epithelial cells can act as effector 
cells, responding to a variety of exogenous and/or en- 
dogenous stimuli by generating and releasing addi- 
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tional mediators of inflammation such as cytokines, 
eicosanoids, and reactive oxygen species, thereby ap- 
pearing to play an important regulatory role in inflam- 
matory and immune processes. Airway epithelial cells 
are known to produce as well as respond to the potent 
inflammatory cytokines interleukin-l(3 (IL-1(3) and tu- 
mor necrosis factor-a (TNF-a). These two cytokines 
have been implicated in the pathogenesis of many acute 
and chronic inflammatory as well as autoimmune dis- 
eases including septic shock, rheumatoid arthritis, and 
inflammatory bowel disease. They are also believed to 
play an important role in the respiratory inflammation 
observed in chronic bronchitis, chronic obstructive pul- 
monary disease, and asthma (14, 15, 40, 45). 

TNF-a and IL-1(3 are produced by a wide variety of 
immune and structural cells including monocytes, 
macrophages, T cells, and endothelial and epithelial 
cells, and they exert many shared biological effects on 
a number of different target cells. They appear to 
activate the same set of transcription factors, primarily 
nuclear factor-ncB (NF-kB) and activator protein-1, 
even though they signal through distinct receptors 
(16). In turn, these two transcription factors govern the 
expression of many genes that are involved in the 
immune responses and inflammatory conditions. 
Genes responding to the effect of TNF-a and IL-1(3 via 
the action of these transcription factors include cyto- 
kines, growth factors, cell adhesion molecules, metal- 
loproteinases, and other inflammatory mediators (3, 
28). 

Airway epithelial cells may participate in the devel- 
opment and maintenance of these respiratory diseases 
through different mechanisms. They may initiate and 
help sustain inflammation through the release of in- 
flammatory mediators and by contributing to cell-spe- 
cific recruitment and activation of granulocytes and 
other leukocytes through differential expression of 
proadhesive signals and/or chemotactic factors and cy- 
tokines (7). In addition, during the course of sustained 
inflammation, the airway epithelium may also be in- 
volved in the remodeling process that is typical of 
obstructive lung diseases through the release of growth 
factors and profibrotic signals (30). 
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GENES INDUCED BY INFLAMMATORY CYTOKINES IN THE AIRWAYS 



The identification of genes in which expression is 
altered in the inflamed airway epithelium can help to 
delineate the molecular mechanisms of the inflamma- 
tory processes observed in respiratory diseases. Our 
recent studies have focused on categorizing the expres- 
sion pattern of transcripts and associating the changes 
to cellular events that ultimately contribute to the 
determination of the respective function of the identi- 
fied genes. cDNA representational difference analysis 
(cDNA-RDA) is a sensitive, PCR-based subtractive en- 
richment procedure that can be used for the isolation of 
genes with altered expression between various tissues 
or cell samples (23, 32). We performed cDNA-RDA to 
identify the genes differentially expressed between 
TNF-a- and IL-l(3-stimulated primary human bron- 
chial epithelial cells (HBECs) and unstimulated con- 
trol cells. 

Here we describe the identification of the transcripts 
more abundantly expressed in HBECs stimulated with 
TNF-a and IL-1(3. Several of these genes have previ- 
ously been reported as being more actively transcribed 
with cytokine stimulation in various cell and tissue 
types (6, 8, 22, 26, 27, 29, 33-35, 37-39, 43). We also 
identified two uncharacterized expressed sequence 
tags (ESTs), one of which showed preferential expres- 
sion in the adult lung along with a sequence not rep- 
resented in the public databases that is predominantly 
expressed in adult and fetal lungs. Identification of the 
genes central to epithelial airway inflammation may 
provide novel targets that may serve to modify the 
inflammatory condition and, ultimately, the progres- 
sion of respiratory disease. 

MATERIALS AND METHODS 

Cell culture. Primary HBECs obtained from Clonetics were 
seeded at a density of 2.5 X 10 5 cells/ml and grown to 70-80% 
confluence in the medium supplied (Clonetics). Actively 
growing cells were simultaneously stimulated with TNF-a 
(100 ng/ml) and IL-1(3 (100 ng/ml) for 4 h (9). Control cells 
were incubated in medium for the same length of time, and 
then both cell populations were harvested. 

RNA isolation and cDNA synthesis. Total RNA was ex- 
tracted with TRIzol reagent (GIBCO BRL, Life Technologies) 
as described by the manufacturer. mRNA was purified with 
the MessageMaker reagent assembly kit (GIBCO BRL) with 
conditions recommended by the manufacturer. For cDNA- 
RDA, 300 ng of mRNA were used to synthesize double-strand 
cDNA with the Superscript choice system (GIBCO BRL). For 
RT-PCR analysis, first-strand cDNA was synthesized from 5 
jig of total RNA that had been DNase treated with RNase- 
free DNase I (Roche Diagnostics). 

Generation of subtracted libraries by cDNA-RDA. cDNA- 
RDA was performed essentially as described by O'Neill and 
Sinclair (32). Briefly, double-strand cDNA was extracted once 
with an equal volume of phenol-chloroform-isoamyl alcohol 
and precipitated with 0.5 volume of 7.5 M ammonium acetate 
and 2.5 volumes of ethanol. Molecular biology-grade chemi- 
cals were purchased from Sigma, and the cDNA precipitation 
conditions used were as recommended by the manufacturer 
of the cDNA synthesis kit (GIBCO BRL). After centrifuga- 
tion, the resultant pellet was washed with 70% ethanol and 
resuspended in 16 jil of sterile water. Double-strand cDNA 



was digested with Dpnll (New England Biolabs) and ligated 
to R-adaptors (a 24-mer annealed to a 12-mer). Amplicons for 
both "tester" (cDNA prepared from stimulated HBECs) and 
"driver" (cDNA prepared from unstimulated control HBECs) 
were generated with Expand Long Template PCR system 
and Expand Buffer 1 (Roche Diagnostics). Typically, five 
100-jxl PCRs for each tester and 10 PCRs for each driver were 
performed. Dpnll digestion was used to remove the R-adap- 
tors from both driver and tester amplicons followed by liga- 
tion of J-adaptors to the tester. Subtractive hybridizations 
were performed in 5-|xl reactions at 67°C for 20 h in a 
GeneAmp 2400 (PerkinElmer). To generate difference prod- 
uct one (DPI), 250 ng of tester cDNA were mixed with 25 |xg 
of driver cDNA at a ratio of 1:100. DPI cDNA was digested 
with Dpnll to remove J-adaptors before ligation of N-adap- 
tors. To generate difference product two (DP2) from DPI, 
31.25 ng of tester were mixed with 25 jig of driver cDNA at a 
ratio of 800:1. Digested and excess adaptors were removed by 
washing the cDNA on Microcon 30 filters (Amicon) as recom- 
mended by O'Neill and Sinclair (32). 

Cloning of difference products. DP2 was amplified with the 
primer 5 '-CUACUACUACUACUAAGGCAACTGTGCTATC- 
CGAGGGAA-3' (N-Bgl-24dU) and cloned into the pAMPIO 
vector with the ligation-independent cloning protocol pro- 
vided with the CloneAMP kit (GIBCO BRL). Plasmids were 
purified from 3-ml cultures of the white colonies with the 
Wizard Plus Minipreps DNA purification kit (Promega). 

DNA sequencing and analysis. Cycle sequencing of plas- 
mid DNA was performed on an automated ABI310 or ABI377 
sequencer (PerkinElmer) with M13 forward and reverse 
primers according to the manufacturer's instructions. Se- 
quence similarity searches were performed with the BLAST 
algorithm (2), and sequence alignments were generated with 
the GCG software package (Wisconsin Package version 9.1). 

Northern and "virtual" Northern blots. For Northern blot 
analyses, a master mix of mRNA or total RNA from unstimu- 
lated and cytokine-stimulated HBECs was combined with 
RNA loading buffer (Sigma), and aliquots were loaded for the 
generation of eight 1-jxg mRNA blots and eight 20-|xg total 
RNA blots. RNA was electrophoresed in 1% agarose formal- 
dehyde denaturing gels for 2 h. RNA was transferred to 
Hybond-N membranes (Amersham) by capillary transfer 
with 20 X saline-sodium citrate (SSC), and ultraviolet cross- 
linked to the filter. Virtual Northern blots (= Southern blots 
of cDNA) were prepared with the same unsubtracted cDNA 
representations that were used as the tester or driver during 
the enrichment procedures. Five hundred nanograms of am- 
plicon cDNA were electrophoresed in 1.5% agarose gels and 
blotted by capillary transfer. 

Hybridization analyses. Clone inserts for the hybridiza- 
tion analysis were generated by PCR amplification from 
plasmid DNA with the N-Bgl-24 primers (32). Probes for 
DP2 cDNA and amplified plasmid inserts were synthesized 
by random priming of 25 ng of DNA in the presence of 
[a- 32 P]dCTP (Amersham) with the Rediprime random- 
primer labeling system (Amersham). All hybridizations 
were performed overnight in ExpressHyb hybridization 
solution (Clontech) at 65°C. After hybridization, the Northern 
and RNA Master blots (Clontech) were washed in 2X SSC- 
0.05% SDS for 20 min at room temperature followed by two 
20-min washes in 0.1X SSC-0.1% SDS at 50°C. For Southern 
blots and GeneFilters cDNA arrays (Research Genetics), 
three consecutive 15-min washes were performed at 65 °C in 
solutions containing 0.5X SSC-0.1% SDS followed by 0.2X 
SSC-0.1% SDS and finally 0.1X SSC-0.1% SDS. The filters 
were exposed to phosphorimager plates for between 2 h and 
5 days and were visualized by a STORM 840 phosphorimager 
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(Molecular Dynamics). Representative HBEC Northern and 
Southern blots were hybridized with a glyceraldehyde-3- 
phosphate dehydrogenase (GAPDH) probe (Ambion) for nor- 
malization of transcript intensities. Individual blots were 
stripped and reprobed no more than 3 times to ensure that no 
significant reduction of signal intensity occurred. Multiple 
tissue Northern blots were hybridized with a GAPDH probe 
(Ambion) to control for loading. GeneFilters hybridization 
results were analyzed with the Pathways software (Research 
Genetics). 

Semiquantitative RT-PCR analysis. Semiquantitative RT- 
PCR analysis for GAPDH, IL-8, the Epstein-Barr virus 
(EBV)-induced gene 3 (EBI3), EXER107, and EXER102 was 
performed. Briefly, 100 ng of single-strand cDNA were used 
to seed a 50-ul PCR mixture and were subjected to 22-36 
cycles of PCR (95°C for 20 s, 58-60°C for 20 s, and 72°C for 
20 s) in a PerkinElmer DNA thermal cycler. After the desig- 
nated number of cycles, 7 jjlI of the reaction mixture were 
removed before the samples were returned for further cy- 
cling. The PCR products were visualized on 1.5% agarose gels 
stained with 0.5 |xg/ml of ethidium bromide. Scanned gel 
images were imported into ImageQuant 5.0 (Molecular Dy- 
namics), and the volume density of the amplified products 
was calculated. 



RESULTS 

cDNA-RDA using unstimulated and cytokine-stimu- 
lated HBEC mRNA. Cultured primary HBECs were 
used to mimic in vivo inflammation and to identify 
cytokine-responsive genes with potential relevance to 
inflammatory airway diseases. Initially, we deter- 
mined the optimal stimulation period for the HBECs 
cultured with TNF-a and IL-1(3. After stimulation, 
RNA was harvested after 1, 2, 4, 8, and 24 h of incu- 
bation. Semiquantitative RT-PCR relative to a GAPDH 
control was performed with several genes, including 
those for IL-8, intercellular adhesion molecule- 1 
(ICAM-1), IL-6, and TNF-a, known to be important in 
airway inflammation and induced by TNF-a and IL-1(3 
(27, 34, 43). Optimal stimulation was determined as 
being 4 h after the initial addition of the two cytokines 
(results not shown). cDNA-RDA was then performed 
with unstimulated and cytokine-stimulated HBEC 
mRNA samples. To isolate genes more actively tran- 
scribed in cytokine-stimulated HBECs, we used cDNA 
from cytokine-stimulated HBECs as the tester and 
cDNA from unstimulated HBECs as the driver. Two 
rounds of enrichment were used to generate DP2 be- 
cause these conditions have previously been shown to 
be sufficient for removing unwanted background and to 
enrich for differentially expressed fragments (32). To 
determine whether the subtraction experiment had 
been successful in enriching for sequences expressed at 
a higher abundance in the tester population, we hy- 
bridized the DP2 population to high-density cDNA 
arrays. Each of the two cDNA filters (GeneFilters 
GF200 and GF202) contains -5,000 arrayed se- 
quences. For GF200 and GF202, -50 and 10%, respec- 
tively, of these sequences were highly similar to or 
represent genes of known function, the remaining se- 
quences on the filters represent ESTs with unknown 
function. To determine which genes were potentially 



represented in DP2, hybridization results were ana- 
lyzed with the Pathways software (Research Genetics). 
Table 1 summarizes these results and demonstrates 
that many of the known genes detected on these filters 
have previously been shown to be upregulated by 
TNF-a and/or IL-1(3 in a variety of other cell and tissue 
types. 

Sequence analysis of subtraction products. The cDNA 
species present in DP2 that had been enriched for 
genes that were upregulated by cytokine stimulation 
were cloned into the pAMPIO vector and further char- 
acterized by DNA sequence analysis. To determine the 
redundancy of the library, inserts from 72 clones were 
amplified and digested with the restriction enzymes 
Rsal and Alul and analyzed on a 1.5% agarose gel. 
After inspection of the digestion products, 48 clones 
were chosen for sequencing. Comparison of these se- 
quences with those in the public databases resulted in 
the identification of 18 known and novel genes (results 
not shown). To determine the redundancy in the li- 
brary, 384 clones from this library were arrayed on 
filters, and probes representing these 18 genes were 
individually hybridized to the arrays. Results from this 
analysis indicated that only 88 clones from this 384- 
clone library (-23% of the library) gave hybridization 
signals. These data indicated that the DP2 cDNA pop- 
ulation potentially contained gene fragments repre- 
senting many more genes. We therefore generated and 
sequenced a 1,152-clone cDNA library. Sequence infor- 
mation was collected, vector sequences were removed, 
and the remaining sequences were assembled into 
clusters. Sequences representing the individual clus- 
ters were compared with those present in the public 
databases. Sixty-six different sequences were repre- 
sented by the 1,152 clones. Of these, 57 were identical 
to known genes, 1 sequence (KIAA0058) identified a 
protein of unknown function, 6 identified human ESTs 
(EST1-6), and 2 sequences did not show any similarity 
to any known sequences (EXER107 and UNK2). Table 
2 summarizes the results of this analysis. Twelve 
clones contained the sequence from the R-Bgl-24 adap- 
tor, the sequence used in the generation of the original 
representations. These clones most probably arose 
from an inefficient digestion with the Dpnll restriction 
enzyme, and hence the R-adaptors were not totally 
removed, resulting in cDNA fragments that would be- 
have effectively as tester unique fragments in the sub- 
tractive hybridization step. Clones identified as con- 
taining the R-adaptor sequence were not analyzed 
further. 

Expression analysis of individual transcripts. Be- 
cause methods based on subtractive enrichment proce- 
dures have been shown to yield false positives, the 
expression pattern of individual clones was further 
verified by virtual Northern and Northern blotting and 
semiquantitative RT-PCR analysis. For virtual North- 
ern and Northern analyses, cDNA inserts were isolated 
from the sequenced clones and used as hybridization 
probes on blots containing unstimulated and stimu- 
lated cDNA and mRNA (Figs. 1 and 2). Not all se- 
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GENES INDUCED BY INFLAMMATORY CYTOKINES IN THE AIRWAYS 
Summary of genes identified by hybridization ofDP2 to GeneFilters GF200 and GF202 cDNA arrays 











Proteases and antiproteases 






225.6 


1 


F 


3 


2 


Hs.75557 


22040 


T72581 










9 


Hs.22239 


30850 


R42600 


230.4 


1 


G 


7 


10 


Hs.79572 


264117 


N20475 










26 


Hs.84898 


898035 


AA598950 


263.7 


1 


E 


2 


11 


Hs.1417 


240766 


H80215 


243.8 


1 


G 


2 


22 


Hs.75621 


123561 


R00822 


258.3 


2 


G 


4 


2 


Hs.14804 


109316 


T80924 










Adhesion molecules 






235.7 


1 


F 


12 


5 


Hs.79133 


40751 


R56219 


266.0 


1 


F 


12 


6 


Hs.76206 


69672 


T53626 


248.5 


1 


G 


6 


7 


Hs.82004 


251019 


H97778 


224.7 


1 


F 


11 


27 


Hs.58464 


774754 


AA442092 










26 


Hs. 82848 


149910 


H00756 


315.9 


1 


A 


2 


27 


Hs.89546 


186132 


H39560 


236.6 


1 


A 


3 


24 


Hs.832 


340644 


W56709 


248.3 


1 


H 


8 


29 


Hs.851 


773332 


AA425451 


253.3 


1 


G 


2 


29 


Hs.78146 


130541 


R22412 


227.5 


1 


D 


10 


11 


Hs.20313 


180298 


R35257 










Transcription factor 






489.1 


1 


A 


7 


25 


Hs.1665 


23804 


R38383 


221.8 


1 


F 


6 


19 


Hs.765 


126368 


R06446 










Cytokines 








167.0 


1 


D 


2 


7 


Hs. 126256 


324655 


W47101 


229.6 


1 


A 


2 


22 


Hs.340 


768561 


AA425102 


430.4 


1 


A 


2 


16 


Hs.69744 


840753 


AA486072 










Membrane pro 








231.8 


1 


A 


3 


26 


Hs.46459 


755054 


AA482637 


462.7 


1 


E 


3 


12 


Hs.23582 


809938 


AA454810 


225.7 


1 


A 


7 


24 


Hs.82685 


813552 


AA455448 


274.5 


1 


H 


12 


2 


Hs.25856 


739183 


AA421296 


260.8 


1 


B 


12 


3 


Hs.46437 


842860 


AA486393 










26 


Hs. 106671 




R56100 










Apoptosis 








238.5 


1 


D 


10 


24 


Hs.38768 


814478 


AA459263 










Growth factors 






223.8 


1 


B 


2 


9 


Hs.73793 


34778 


R19956 










Unassigned 






5757.2 


1 


D 


7 


14 


Hs.75522 


810444 


AA457114 


286.1 


1 


G 


3 


15 


Hs.25590 


547247 


AA085318 



quences were tested for differential expression because 
many genes have been previously reported as having 
induced expression with TNF-a and IL-ip. We tested 
31 sequences by virtual Northern analysis for the dif- 



ferential abundance between tester and driver cDNAs. 
We detected an increased hybridization intensity in 
cDNA generated from cytokine-stimulated HBECs for 
23 of these sequences (Fig. 1), 3 of the sequences did 
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Table 1. — Continued 



Upregulation 
Reported, 

Title Ref. No. 




Chromosome 


STS Name 


Proteases and antiproteases 






Matrix metalloproteinase-2 (gelatinase A; collagenase type IV) 6*t 


MMP2 


16ql3 


WI-7659 


Homo sapiens mRNA for putative MT4-MMP protein (GF202) 










CTSD 


llpl5.5 




Cathepsin B$ ? 26* 


CTSB 


8p22 


SHGC-9711 




TIMP1 






activity, collagenase inhibitor) 










PI 


14q32.1 


SHGC-32639 










(H. sapiens)^. 
















Cadherin 8 Adhesion molecules 


CDH8 




SHGC-15929 


Cadherin 5, VE-cadherin (vascular epithelium) 


CDH5 


16q22.1 


SHGC-13798 


H. sapiens mRNA for E-cadherin 






SHGC-11833 


Catenin (cadherin-associated protein), pi (88 kDa) 


CTNNB1 


3p22-p21.3 




L-selectin (lymphocyte adhesion molecule-1) 






WI-7045 


E-selectin (endothelial adhesion molecule-1) 38 t 


SELE 


Iq23-q25 


WI-9125 




ITGB8 


7 


sWSS2601 


Integrin a E (antigen CD103 human mucosal lymphocyte 


ITGAE 


17 




antigen-1; a-polypeptide) 








Platelet endothelial cell adhesion molecule (CD31 antigen) 


PECAM1 


17q23 


WI-7753 


Human cell adhesion kinase-p (CAKp) mRNA, complete cDNA 




8 


WI-12748 


sequence (cds) 








Transcription factor 








Z' fi t ' h 1 t Zf 36 ' 8* 
PATAKWl^ 10 eU1 t 0m i ° f °?°^ S ? P " • 1 ^' In 'i US ^ *s 


ZFP36 






in ing pro ein- (g o in ranscrip ion ac or- ) 




Xpll.23 




Cytokines 








Interleukin-ip (GF202*)$ 




2ql3-q21 


CHLC.UTR 00699 X04500 




SCYA2 


17qll.2-ql2 


WI-7605 


S ™ a ho i mo"ogous°to mou.se^ig'je) 1005 '* 6 cllemoatt ' ract ' aIlt P rotein 22 








Small inducible cytokine A5 (RANTES) (GF202) 35*t 


SCYA5 


17qll.2-ql2 


CHLC.UTR_03943_M21121 


Membrane proteins 








Human putative transmembrane receptor IL-1 receptor-related 








protein mRNA, complete cds [interleukin-18 receptor 1 








(IL-18R1)] 








Membrane component, chromosome 1, surface marker 1 (40 -kDa 


M1S1 


Ip32-ql2 


W 1-7569 


glycoprotein, identified by monoclonal antibody GA733) 








CD47 antigen (Rh-related antigen, integrin-associated signal 


CD47 


3ql3.1-ql3.2 


WI-8002 


transducer) 








CD68 antigen 


CD68 






Cytokine receptor family II, member 4 


CRFB4 


21q22.1 


CHLC.UTR_04215_Z17227 


ESTs, highly similar to thymic epithelial cell surface antigen 








(Mus musculus) (GF202) 








Apoptosis 








Human Bcl-2 related (Bfl-1) mRNA, complete cds$ 29*t 






SHGC-35552 


Growth factors 








Vascular endothelial growth factor 37* 


VEGF 


6pl2 


SHGC-9780 


Unassigned 








B94 proteinf 39* 




14 




Human stanniocalcin precursor (STC) mRNA, complete cds$ 




8 


STS_U25997 



DP2, difference product 2; MT4-MMP, membrane type 4 matrix metalloproteinase; EST, expressed sequence tag; RANTES, regulated on 
activation normal T cell expressed and secreted. References are for previous studies that reported upregulation of the transcript or protein 
in airway epithelial cells or other cell systems by tumor necrosis factor-a (*) or interleukin (IL)-ip (t). $ Fragments from these genes were 
also identified in the sequenced library cloned from DP2. 



not show appreciably different hybridization intensi- 
ties between the 2 lanes, and 5 sequences did not 
generate a hybridization signal (Table 2). Nineteen 
sequences were analyzed by Northern blotting, and 



differential expression was shown for 10 known genes, 
2 ESTs (EXER102 and EST3), and a novel sequence 
(EXER107; Fig. 2) compared with a GAPDH control. 
Six sequences did not show any obvious visible expres- 
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Summary of genes identified by sequencing of DP2 library 



Differential Differential 
Expression Expression 



Gene Product 


Virtual 
Northern 


Northern 


Gene Product 


Virtual 

Northern Northern 


Membrane proteins 






; 

Housekeeping 














H. sapiens syndecan 4 (amphiglycan, 


Yes 


"ies 


H. sapiens [3-actin (ACTB) 




ryudocan) (SDC4) 










Human cytokine receptor (EBI3)* 


ND 


Yes 


fJ-Tubulin 




Annexin VIII (ANX8) 


Yes 




P-Tubulin-like 




Cytokines/chemokines/secreted proteins 






Lactate dehydrogenase-A (LDH-A) 




ENA-78 


Yes 


Yes 


H. sapiens aldehyde dehydrogenase 6 




Human follistatin 


Yes 


Yes 


(ALDH6) 




IL-ipt 


Yes 




Human acidic ribosomal 




Tumor necrosis factor (TNF-a) 






phosphoprotein PO mRNA 




Granulocyte-macrophage colony-stimulating 










factor (GM-CSF) 






Keratins 




Granulocyte colony-stimulating factor (G-CSF) 






Human type II keratin K5 




Proteases/antiproteases 






H. sapiens keratin 6 




MMP1 


Yes 




Human mesothelial keratin K7 




MMP3 




Yes 


Keratin K19 




MMP9 










MMP10 






EXER102 


ND Yes 


Urokinase-type plasminogen activator 


Yes 


Yes 




Yes Yes 


Tissue-type plasminogen activator (tPA) 


Yes 


NS 






H. sapiens mRNA for trypsinogen IV a-form 








ND 


H. sapiens ai-antichymotrypsint 


Yes 


NS 


Protein of unknown function 




Cathepsin Bf 






Human mRNA for KIAA0058 


Yes NS 


Apoptosis 










Human Bcl-2 related (Bfl-1), BCL2A1, GRSf 








V 


Human lamin B2 (LAMB2) 


Yes 


NS 


Clones containing R-Bgl-24 sequence 




Kinases 






Human membrane glycoprotein 4F2 




H. sapiens protein kinase, mitogen-activated 










kinase 3 






protein 2 (CRABP2) 




Human mRNA for phosvitin/casein kinase II (3 


Yes 


NS 


Transcriptional intermediary factor-l[3 




Unassigned 






zinc finger protein, nuclear 




B94f 


Yes 


Yes 


corepressor KAP-1 




Stanniocalcin (STCl)t 


Yes 




H. sapiens rabkinesin6 




H. sapiens B-factor, properdin (BF), 










complement B 






inhibitor (kop) 




H. sapiens X (inactive )-specific transcript 


ND 














HumanTsUgene from interferon 




Human thyroid hormone binding protein (p55) 


ND 




inducible gene 




Platelet-type phosphofructokinase 






H. sapiens fibulin 1 (FBLN1) 




Human transglutaminase mRNA 






H. sapiens phosphatidylinositol 




Stress/antioxidant proteins 






synthase (PIS) 




Human 90-kDa heat shock protein 


Yes 


Yes 


EST1 




Mn superoxide dismutase 


Yes 




EST6 




Human glutathione S-transferase (Plb) 






UNK2 




Signal transduction 










Human epithelial cell marker protein 1 


NS 








(HMel) 14-3-3 a-protein 










H. sapiens moesin (MSN) 


NS 









ICAM-1, intercellular adhesion molecule-1; ENA-78, epithelial cell-derived neutrophil-activating protein-78; ND, not detected; NS, not 
shown. * Genes tested by RT-PCR. t Genes also identified by hybridization of DP2 to the arranged cDNAs. 



sion difference between unstimulated and cytokine- 
stimulated HBEC mRNAs (Table 2). All of the genes 
that showed differential expression by Northern blot- 
ting had an increased abundance in tester compared 
with driver cDNA when analyzed by virtual Northern 
blotting except for EBI3 and MMP3, which were not 
detected or tested, respectively, on virtual Northern 
blots (Table 2). We were not able to detect significant 
differential expression for five transcripts that had an 
increased abundance in tester compared with driver 
cDNA when analyzed by Northern blotting. This sim- 



ply may be the consequence of transcript levels not 
having been increased sufficiently for visible detection 
(Table 2). The data for the transcripts that were de- 
tected by Northern blotting indicated that — 70% of the 
genes identified by the cDNA-RDA experiments were 
indeed differentially expressed. To confirm differential 
expression for several genes in which expression was 
only weakly detected by Northern blotting, we per- 
formed semiquantitative RT-PCR relative to GAPDH 
transcript levels. These data indicated that the genes 
represented by the clones for EXER107, EBI3, and 
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Mitogen activated 
kinase 3 

Phosvitiit casein 
II beta 




Tissue plasminogen 
activator 
Urokraase-type 
plasminogen activator 



9(lkDa Heat 
Shock Protein 




GA.P.DH 



Fig. 1. Virtual Northern analysis of individual clones. Inserts from individual clones isolated from difference 
product 2 (DP2) were used as hybridization probes on virtual Northern blots containing cDNA from unstimulated 
( — ) and tumor necrosis factor (TNF)-a- and interleukin (IL)-lpi-stimulated ( + ) human bronchial epithelial cells 
(HBECs). Five hundred nanograms of tester and driver representations were loaded in each lane. Loading of 
Northern blots was corrected for by hybridizing a glyceraldehyde-3-phosphate dehydrogenase (GAPDH) probe to 
representative blots. ICAM-1, intercellular adhesion molecule-1; MMP1, matrix metalloproteinase 1; Mn SOD, 
manganese superoxide dismutase; EST3 and EST4, expressed sequence tags 3 and 4, respectively. See text for 



EXER102 were expressed more abundantly in cyto- tial expression; however, our data clearly indicate the 

kine-stimulated HBECs (Fig. 3). need to perform Northern analysis or quantitative 

Several reports (17, 24, 46) have used hybridization RT-PCR for confirmation of differential expression. Of 

to cDNA representations as an indication for differen- 14 genes that clearly showed differential hybridization 



90kDa Heat 
Shock Protein 



(Bcl-2 related) 



Mitogen activated 
kinase 3 



Urokinase type 
plasminogen activator 



3 



Fig. 2. Northern analysis of individual clones. Inserts from individual clones derived from DP2 were used as 
hybridization probes on Northern blots containing 10 u.g of total RNA from unstimulated and TNFa- and 
IL-ip-stimulated HBECs. Loading of virtual Northern blots was corrected for by hybridizing a GAPDH probe to 
representative blots. ENA-78, epithelial cell-derived neutrophil-activating protein-78. See text for discussion. 
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EXER187 (30-34) 





|i: 







+ 



Fig. 3. Semiquantitative RT-PCR was performed 
relative to GAPDH amplification. Left: 100 ng of 
single-strand cDNA from unstimulated (unstim) 
and TNF-a- and IL-ip-stimulated (stim) HBECs 
were used as templates in both reactions. Seven 
microliters of amplification product were analyzed 
on a 1.5'/i agarose gel. Nos. in parentheses, cycle no. 
Right: individual amplification products were ana- 
lyzed for integrated pixel intensity (jy-axis) and plot- 
ted against cycle number (x-axis) with ImageQuant 
5.0 software. As a control for HBEC cytokine stim- 
ulation, an IL-8 amplification was performed. EBI3, 
Epstein-Barr virus-induced gene 3. See text for dis- 





GAPDH (22-26) 



^I$PP : iiiiiii 



by virtual Northern analysis, only 9 showed distinct 
visible differential expression, whereas 5 genes did not 
when tested on Northern blots. 

Tissue distribution of EST3, EXER102, and 
EXER107 expression. Novel genes represented by the 
clones for EST3, EXER102, and EXER107 were fur- 
ther analyzed with multiple tissue Northern blots 
and RNA master blots. The probe for EXER102 de- 
tected transcripts of ~3 and 6 kb in the majority of 
mRNAs from adult and fetal tissues (Fig. 4). The 



most abundant expression was found in adult and 
fetal kidneys and adult lung. Analysis of the expression 
of EXER102 with the RNA Master blot indicated that 
the comparatively highest levels of expression were in 
the trachea, lung, small intestine, kidney, pituitary 
gland, stomach, thyroid, and fetal lung (results not 
shown). The probe for EXER107 detected a single tran- 
script of —3.5 kb expressed at varying degrees in all 
adult and fetal tissues; however, a significantly higher 
abundance of the message was found in adult and fetal 
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Fig. 4. Northern blot analysis for EXER102 
and EXER107. Inserts from representational 
difference analysis (RDA) clones were hybrid- 
ized to multiple tissue Northern blots. To 
control for differences in loading between 
lanes, the same filters were hybridized with a 
GAPDH control. mRNA from indicated tis- 
sues was loaded in each lane. 



GAPDH 



lungs (Fig. 4). Analysis of the expression of EXER107 
with the RNA Master blot indicated that the compar- 
atively highest levels of expression were found in the 
trachea, adult lung, kidney, stomach, and thyroid (re- 
sults not shown). Hybridization of the EST3 probe to 
multiple tissue Northern blots did not result in the 
detection of a transcript in any tissue tested; however, 
hybridization to a RNA Master blot indicated that the 
gene represented by EST3 was expressed most abun- 
dantly in the pancreas, mammary gland, thyroid, blad- 
der, and trachea (results not shown). 

DISCUSSION 

We performed cDNA-RDA to identify and isolate 
genes in primary HBECs in which expression is 
induced by the inflammatory cytokines TNF-a and 
IL-1(3. We used a mixture of these two cytokines in 
our experiments because their coexpression and 
their synergistic effects are important for disease 
pathology in inflammatory conditions. As expected, 
we found an induced expression of a number of 
known inflammation-related genes, including IL-1(3 
and TNF-a themselves; several genes involved in 
inflammatory cell trafficking, such as IL-8, regulated 
on activation normal T cell expressed and secreted 
(RANTES), and epithelial cell-derived neutrophil- 



activating protein-78 (27, 34, 43); and genes involved 
in the activation of inflammatory cells, such as gran- 
ulocyte-macrophage colony-stimulating factor and 
granulocyte colony-stimulating factor, that confirms 
the validity of our approach. In addition, the stimu- 
lated airway epithelial cells showed induced expres- 
sion of a number of protective genes including those 
for antiproteases, such as the tissue inhibitors of 
metalloproteinases (5), and antioxidants, such as 
manganese superoxide dismutase (47). TNF-a is also 
known for its ability to induce apoptosis in a variety 
of cell types through its activation of its cognate 
receptor subtype I (18). One mechanism by which the 
airway epithelial cells may protect themselves against 
TNF-a-induced apoptosis is by expressing Bfl-l/Al, a 
protein that is able to block cellular apoptosis. This 
prosurvival Bcl-2 homolog is a direct transcriptional 
target of NF-kB, hence providing evidence supporting 
the protective role of NF-kB in the immune system 
(51). Another gene involved in apoptosis regulation 
that we identified as being upregulated by TNF-a and 
IL-1(3 is the gene for the 90-kDa heat shock protein 
(HSP90). HSP90 functions as a molecular chaperone 
that aids in the folding of nuclear hormone receptors 
and protein kinases. A study (19) has shown that an 
excess of HSP90 is associated with increased apoptosis 
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and, conversely, that reduced HSP90 levels correlate 
with protection against apoptosis in TNF-a- and cyclo- 
heximide-treated cells (19). It therefore seems possible 
that in bronchial epithelial cells exposed to inflamma- 
tory conditions, a similar fine balance between survival 
and apoptosis may exist. In vivo, when this balance is 
shifted in favor of apoptosis, epithelial cell shedding 
may occur as seen in the asthmatic condition (25). 

We also identified two kinases involved in intra- 
cellular signaling as having upregulated transcrip- 
tion on cytokine stimulation. Previous work (48) has 
shown that mitogen-activated protein (MAP) kinase 3 
(MKK3) is a critical component of a TNF-a-stimulated 
signaling pathway that causes increased expression of 
inflammatory cytokines. Targeted disruption of the 
MKK3 gene resulted in a selective defect in the re- 
sponse of fibroblasts to TNF-a, including reduced p38 
MAP kinase activation and cytokine expression (48). 
However, its upregulation in inflamed airway epithe- 
lial cells has not been previously reported. Phosvitin/ 
casein kinase II (CKII) was also shown to be present in 
DP2 (Table 2 ). This kinase has previously been shown 
to immunoprecipitate as part of the NF-kB cytoplasmic 
complex, and phosphorylation of NF-kB by CKII may 
control translocation of the activated complex or mod- 
ulate its ability to bind DNA (4, 44). 

The cell surface proteoglycan syndecan-4 was iden- 
tified as being upregulated in HBECs by cytokine stim- 
ulation. This molecule has the ability to regulate cell 
adhesion via its glycosaminoglycan chains and the 
discrete domains of its core protein; it is also able to 
bind and activate protein kinase C-a (21). Previous 
work in vitro (50) has shown that TNF-a is able to 
induce endothelial cell syndecan-4 expression by both 
increasing gene expression in a NF-HcB-dependent 
manner and prolonging mRNA half-life. In vivo studies 
have shown that levels are also increased after various 
forms of tissue injury including skin wounds, vascular 
wall injury, or myocardial infarction (50). Both RT- 
PCR and Northern analysis (results not shown) indi- 
cated that the transcript levels for another cell surface 
molecule, the Epstein-Barr virus (EBV)-induced gene 3 
(EBI3), were also increased in cytokine-stimulated 
HBECs. EBI3 encodes a 34-kDa glycoprotein, a mem- 
ber of the hematopoietin receptor family related to both 
the p40 subunit of IL-12 and the ciliary neurotrophic 
factor receptor. The protein lacks a membrane- anchor- 
ing motif and is secreted but is also present on the 
plasma membrane of EBV-transformed B lymphocytes 
and transfected cells (12). Latent membrane protein 1 
of the EBV has been shown to induce the expression of 
EBI3 along with ICAM-1, lymphocyte function-associ- 
ated antigen-3, and CD40 in EBV-negative Burkitt 
lymphoma cells and in epithelial cells via NF-kB acti- 
vation. (13). Consistently, our data indicate that this 
gene can also be induced in HBECs by the two cyto- 
kines TNF-a and IL-1(3, suggesting a possible role for 
this molecule in respiratory inflammation. 

Follistatin is a secreted protein that is able to bind 
and neutralize the actions of many members of the 
transforming growth factor-(3 family of proteins. Fol- 



listatin is known to bind activin, a member of this 
family, and both are upregulated by mediators of in- 
flammation (36). Our data are, however, the first to 
suggest that follistatin plays a possible role in airway 
inflammation. The urokinase-type plasminogen activa- 
tor (uPA) and the uPA receptor (uPAR) are thought to 
be critical in the cell invasion processes. The balance 
between uPAR-bound uPA and its inhibitors modu- 
lates a pericellular proteolytic activity able to give 
"stop-and-go" signals to invading cells, possibly in con- 
cert with the matrix metalloproteinase system. Other 
lines of evidence suggest that the uPA-uPAR interac- 
tion elicits chemotaxis, chemoinvasion, and cell multi- 
plication, with uPA also being implicated in the inflam- 
matory disease rheumatoid arthritis (11). Our data 
now suggest that the uPA may also have a role in lung 
inflammation. The tissue-type plasminogen activator 
(tPA) gene was also identified as being regulated by 
TNF-a and IL-1(3. tPA has been shown to have anti- 
inflammatory properties, and a study (42) suggested 
that tPA may suppress neutrophil activation in vivo 
but does not reduce lung neutrophil infiltration. In- 
deed, treatment with tPA has been shown to signifi- 
cantly improve the mortality rates of patients with 
acute myocardial infarction, and the potential for tPA 
treatment in other inflammatory diseases such as 
acute respiratory distress syndrome and rheumatoid 
arthritis, in which neutrophil-mediated injury is likely, 
has been suggested (41). Our findings add further 
support for its use as a possible treatment for certain 
lung inflammatory diseases. 

Several cDNA fragments that contained the R-Bgl-24 
adapter sequence were identified from DP2. These 
cDNA fragments would effectively behave as tester- 
unique fragments in the subtractive hybridization 
steps. The most likely explanation for the presence of 
these clones in DP2 may be due to inefficient digestion 
by the Dpnll restriction enzyme. Several cDNA frag- 
ments that represent genes generally classified as 
"housekeeping" were also identified from DP2 (see Ta- 
ble 2). These may also be artifacts of the experimental 
process not being successfully eliminated during the 
subtraction procedure. However, identification of cDNAs 
for housekeeping genes among subtraction products has 
previously been reported by several groups including 
our own (10, 17, 20); however, further rounds of sub- 
traction may result in the elimination of these genes. 
Alternatively, it is also possible that some of these 
genes are actually induced by cytokine treatment and 
may indicate a change in the metabolic activity of 
cytokine-stimulated HBECs (6). 

Our studies demonstrate that hybridization of DP2 
to arrayed cDNAs allowed rapid identification of the 
genes represented and provided information regarding 
the success of the subtraction experiment. This ap- 
proach identified several genes including those for E- 
selectin, Zfp-36, monocyte chemoattractant protein-1, 
vascular endothelial growth factor, and RANTES (Ta- 
ble 1), all having been previously reported to be regu- 
lated by TNF-a and/or IL-1(3 (22, 35, 37, 38); however, 
these genes were not identified when we sequenced our 
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1,152-clone library. Differing cloning efficiencies of in- 
dividually amplified fragments may possibly explain 
the reason for this discrepancy. Although hybridization 
to gene arrays can be used to rapidly identify some of 
the genes represented in a complex (DP2), it is unlikely 
to be able to identify all of the genes represented until 
all the genes within the genome have been identified 
and subsequently arrayed. Our data therefore support 
the notion that hybridization of DP2 to gene arrays 
combined with cloning and sequencing may provide a 
more comprehensive analysis of the subtraction prod- 
uct. 

A full-length uncharacterized gene, KIAA0058, was 
also identified as being represented in DP2. Northern 
blot analysis has previously shown this gene to be 
expressed in a variety of tissues (31). The mouse ho- 
mologue, designated prtb (proline codon-rich tran- 
script, brain expressed) has also been recently identi- 
fied, and knockout mice have been generated. These 
mice were viable and fertile and did not display any 
obvious abnormalities (49). Our data suggest that 
KIAA0058 may play a role in lung inflammation (Fig. 
1); however, a more detailed analysis of the lung func- 
tion and bronchial hyperresponsiveness of the prtb 
knockout mice may provide further insight into its role. 

Northern analysis of the gene fragments represent- 
ing genes of unknown function, EXER107 and EXER102, 
indicated that in adult tissue they are predominantly 
expressed in the lung (Fig. 3), with EXER107 also 
appearing important for fetal lung development (Fig. 
3). Although we detected a hybridization signal in total 
RNA derived from HBECs for EST3 (Fig. 2), unexpect- 
edly, we were not able to detect any signal in adult or 
fetal tissues. This result may be explicable due to EST3 
having a restricted pattern of expression unique to 
HBECs, whereas EXER102 and EXER107 may be ex- 
pressed in other cells of the airways and hence give 
detectable hybridization intensities in complex tissues. 
Sequences of EST3, EXER102, and EXER107 showed 
no similarity to any genes of known function in the 
public databases, and, therefore, we are not able to 
speculate further on their potential roles. Further 
work needs to involve isolation of full-length sequences 
for these clones, and this may provide further insight 
into their action. A thorough and quantitative mapping 
of the gene expression pattern of our panel of genes in 
various (patho)physiological situations will give indis- 
pensable information concerning inflammatory condi- 
tions of the lung. This information, in turn, may help 
us devise new rational and gene therapeutic ap- 
proaches for the treatment of various respiratory dis- 
eases. 

We thank Michael O'Neill and Andrew Sinclair for providing us 
with their detailed protocol for cDNA representational difference 
analysis from small amounts of starting mRNA and also Mike 
Hubank and David Schatz for providing us with their original pro- 
tocol. We thank Jennifer McDonald and Maureen Milnamow for 
high-throughput sequencing of the difference product 2 library and 
Carol Jones for critically reviewing our manuscript. 
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